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Abstract.  In  the  stratified  ocean,  vertical  motions  arise  from  both  adiabatic 
and  diabatic  mechanisms.  Diapycnal  advection  is  the  vertical  component 
of  flow  across  an  isopycnal  surface  which  occurs  when  mixing  produces 
a divergent  flux  of  buoyancy.  Buoyancy  forcing  of  the  lateral  flow  by 
vortex  stretching  occurs  when  diapycnal  advection  rates  vary  with  depth. 
Microstructure  observations  of  enhanced  turbulence  above  rough  topography 
are  presented.  These  data  allow  us  to  distinguish  the  influence  of  diabatic 
forcing  on  the  circulation.  Estimates  of  diapycnal  advection  are  used  to 
quantify  the  diabatic  flow,  and  the  role  of  diapycnal  advection  as  a mechanism 
of  vortex  stretching  on  the  circulation  is  assessed.  The  divergence  of  diapycnal 
mass  flux  is  found  to  be  a significant  forcing  mechanism  for  the  circulation 
occurring  above  fracture  zone  topography.  Moreover,  both  the  direction 
and  magnitude  of  diapycnal  advection  are  dependent  on  the  efficiency  of 
turbulence  at  generating  a buoyancy  flux.  Depth  variations  of  diapycnal 
advection  suggest  the  mixing  efficiency  is  reduced  in  weakly  stratified  abyssal 
canyons. 


Introduction 

The  sun  heats  the  Earth’s  atmosphere  and  ocean, 
driving  the  winds,  evaporation,  and  rain.  These  act  over 
global  scales  to  input  mechanical  energy  and  variance 
to  the  scalar  fields  of  the  ocean.  Turbulence  acts  at  the 
smallest  scales  of  fluid  motion,  where  variance  is  dissi- 
pated by  molecular  processes.  Stirring  is  the  mechanism 
that  acts  at  the  intermediate  range  of  scales,  advecting 
variance  both  in  physical  and  wavenumber  space. 

The  input  of  variance  at  the  largest  scales  is  related 
to  the  mechanical  and  buoyancy  forcing  that  drives  the 
so-called  “large-scale  circulation”.  However,  circula- 
tions can  be  forced  by  mechanisms  occurring  anywhere 
along  the  variance  cascade.  For  example,  Spall  (1994) 
describes  large-scale  abyssal  recirculations  driven  by  the 
eddy  fluxes  of  topographic  waves.  Examples  of  eddy 
driven  mean  flows  are  described  by  Lozier  (1997)  and 
Alves  and  Colin  de  Verdiere  (1999).  Buoyancy  forcing 
may  also  originate  at  any  scale.  In  particular,  turbulent 
mixing  at  the  smallest  scales  of  motion  can  give  rise  to 


divergent  fluxes  of  heat  and  salt.  If  the  buoyancy  forces 
that  result  from  turbulence  are  significant  relative  to 
other  forcing,  a large-scale  circulation  may  develop  by 
the  direct  influence  of  small-scale  mixing.  Spall  (2000) 
demonstrates  that  mixing  near  large-scale  topography  is 
particularly  effective  at  driving  strong  buoyancy-forced 
circulations. 

Stommel  (1957)  first  proposed  that  the  abyssal  circu- 
lation of  the  ocean  may  be  driven  by  buoyancy  forcing, 
and  a model  proposed  by  Stommel  and  Arons  (1960) 
described  a simple  set  of  dynamics  governing  the  buoy- 
ancy driven  flow.  The  Stommel-Arons  model  is  still 
central  to  our  view  of  the  abyssal  circulation  problem, 
and  the  important  role  of  buoyancy  forcing  in  the  abyss 
is  not  questioned.  However,  the  nature  of  deep  buoy- 
ancy forcing  has  only  recently  come  to  light,  with  ob- 
servations providing  new  insight  into  the  mechanisms 
driving  abyssal  circulations.  Observations  of  turbulence 
across  the  Brazil  Basin  show  that  abyssal  mixing  rates 
are  spatially  nonuniform,  with  enhanced  mixing  occur- 
ring above  regions  of  rough  topography  ( Polzin  et  al. , 
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1997).  Ledwell  et  al.  (2000)  and  St.  Laurent  et  al. 
(2001)  suggest  that  energy  from  internal  tides  sustains 
the  elevated  turbulence  in  these  regions. 

In  the  report  that  follows,  the  implications  for  buoy- 
ancy forcing  by  small-scale  mixing  are  discussed.  Obser- 
vations from  the  abyssal  Brazil  Basin  described  by  St. 
Laurent  et  al.  (2001)  are  summarized  here.  Estimates 
of  cross-isopycnal  motion,  termed  diapycnal  advection, 
are  used  to  quantify  the  diabatic  flow.  The  diapycnal 
advection  serves  as  a primary  measure  of  buoyancy  forc- 
ing- 


w.(2u) 


Diapycnal  advection 

We  wish  to  consider  the  diabatic  flow  in  the  ocean 
interior  driven  by  a divergent  buoyancy  flux.  The  equa- 
tion for  buoyancy  (6  = —(g/po)p')  is 


|^  + u-Vb  + wN2  = -V-F-^,  (1) 

where  u = (u,  v)  is  the  lateral  flow,  V = (d/dx,d/dy) 
is  the  lateral  differential  operator,  and  (F,  Jb)  are  the 
lateral  and  vertical  components  of  buoyancy  flux.  To 
distinguish  diabatic  from  adiabatic  flow,  it  is  useful  to 
consider  the  buoyancy  budget  along  a neutral  surface. 
Along  a neutral  surface  whose  depth  is  described  by 
Zi(x,y,t),  buoyancy  is  subject  to  the  relations 


V&x  = V6  + VziN2  — 0, 

dbj  I da  \t2  n 

at  ~ dt  ^ at JV 


(2) 


where  (•)*  denotes  an  isopycnal  function  of  (x,  y,  t),  and 
N2  = db/dz.  Additionally,  since  buoyancy  flux  across 
an  isopycnal  occurs  only  by  small-scale  turbulence,  all 
fluxes  may  be  expressed  in  terms  of  a turbulent  diffu- 
sivity  kp  such  that 


Fi  = — kpVbj  = -kpVb  - kpVz;N2  = 0, 
Jb  = -kpN2. 


(3) 


Thus,  the  equation  for  buoyancy  on  a neutral  surface 
can  be  expressed  as 


Vzi)  N2  = ~ v • {kpVZiN2) . 

(4) 

■ We  define  the  vertical  advection  quantity  w*  — w - 
dzi/dt  — u • Vz;  as  the  “diapycnal  advection”.  As  dis- 
cussed by  McDougall  (1987,  1995)  and  Pedlosky  (1996), 


Figure  1.  Diagram  showing  the  parameters  of  a 
steady-state  density  layer. 


w,  is  the  vertical  component  of  flow  through  (not  nor- 
mal to)  isopycnals,  and  is  the  vertical  motion  that  is 
caused  by  a divergent  buoyancy  flux.  The  divergent 
buoyancy  flux  is  comprised  of  two  terms;  a standard 
vertical  divergence  term  —dJb/dz,  and  a lateral  diver- 
gence of  vertical  buoyancy  flux  V ■ (Vz,  Jb) . It  is  use- 
ful to  consider  the  ratio  of  these  terms,  with  scaling 
dJb/dz  ~ S~lJb,  and  V • (Vz,  Jb)  ~ £~2SJb  for  charac- 
teristic lateral  and  vertical  scales  (£,  S).  The  ratio  of  the 
vertical  to  lateral  divergence  terms  is  (f/d)2.  In  general, 
we  will  assume  that  the  aspect  ratio  5/£  of  isopycnals 
is  small,  such  that  {£/8)2  » 1,  and  (4)  can  be  approxi- 
mated as 


(5) 

Dynamical  Considerations 

Consider  a large-scale  geostrophic  and  inviscid  flow. 
The  potential-vorticity  equation  for  flow  in  a steady- 
state  density  layer  (Fig.  1)  is 

fiv  - j-u-  Vh+  ^ (w„(zu)  - w„(zi)) , (6) 

where  h(x,  y)  is  the  thickness  of  the  density  layer,  and 
(zu,  zi ) denote  the  upper  and  lower  bounding  isopycnal 
surfaces.  The  above  expression  is  the  density-layered 
version  of  the  Stommel  and  Schott  (1977)  beta-spiral 
equation,  and  is  discussed  by  McDougall  (1988)  and 
Hautala  and  Riser  (1993).  In  (6),  the  advection  of  plan- 
etary vorticity  is  balanced  by  vortex  stretching.  Vortex 
stretching  may  occur  adiabatically  by  flow  along  sloping 
isopycnals,  or  diabatically  by  the  divergence  of  diapyc- 
nal advection.  Diabatic  vortex  stretching  is  related  to 
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the  divergence  of  buoyancy  flux  by  (5).  During  mixing 
events,  vertical  exchanges  of  buoyancy  are  accompanied 
by  dissipation  of  turbulent  kinetic  energy  (TKE),  and 
the  energy  budget  for  an  ensemble  of  turbulent  events 
is 

(1  — Rf  )Jb  + Rfe  = 0,  (7) 

where  e is  the  dissipation  rate,  and  Rf  is  the  efficiency 
of  the  mixing.  The  TKE  budget  expressed  in  (7)  as- 
sumes steady-state  and  homogeneous  statistics  of  the 
Reynolds  stress  tensor.  This  equation  is  discussed  by 
Osborn  (1980),  and  an  expression  for  the  eddy  diffusiv- 
ity  is  readily  derived  by  substituting  Jj  = —kpN2  into 
(7), 


where  the  lower-case  variables  are  now  dimensionless. 
The  ratio  of  the  diabatic  to  adiabatic  stretching  is 

diabatic  stretching  _ E 
adiabatic  stretching  U2 ' 

The  relative  importance  of  the  diabatic-forcing  term 
will  be  set  by  the  energy  level  of  the  process  controlling 
the  vertical  mixing.  In  the  case  of  dissipating  internal- 
tide  energy,  the  energy  flux  scales  as  U2idc.  In  the  abyss 
where  geostrophic  flow  is  weak,  typical  geostrophic  ve- 
locities are  U ~ 5 mm  s-1,  while  tidal  velocities  are 
generally  Ut ide  ~ 20  mm  s_1.  Thus  a ratio  of  0(1)  for 
the  diabatic  to  adiabatic  vorticity  forcing  is  easy  to  jus- 
tify, demonstrating  the  strong  dynamical  link  between 
mixing  and  circulation  in  the  abyss. 


The  fraction  T = Rf  (1  — Rf)~x  (equivalently  T = 
—Jb/e),  the  ratio  of  buoyancy  flux  to  energy  dissipated, 
is  typically  taken  to  be  <20%  for  stratified  turbulence. 
With  use  of  (5)  and  (7),  the  diapycnal  advection  can  be 
expressed  as 


w. 


dz  ’ 


(9) 


and  the  potential  vorticity  equation  (6)  can  be  ex- 
pressed as 


0v  = £u  - Vh  + — 


(10) 


where  we  have  taken  (w*(zu)  — w,(zi ))  jh  = dwm/dz. 

As  a clear  way  of  demonstrating  the  importance  of 
the  diabatic  forcing  term  in  (10),  it  is  useful  to  consider 
the  following  scaling  relations.  Using  standard  nota- 
tion, we  scale  (x,  z)  ~ (L,  H),  u ~ U,  /3  ~ U/L2  and  de- 
fine the  deformation  radius  as  Rd  = NH/f.  The  scaling 
for  the  term  V/i  is  taken  from  the  geostrophic  relation 
fU  r>j  g'Vh,  where  the  reduced  gravity  is  g'  ~ N2H. 
Additionally,  we  scale  the  energy  dissipation  term  as 
e ~ fE,  where  E is  the  scale  of  the  kinetic  energy  as- 
sociated with  the  forcing  for  the  vertical  mixing.  This 
scaling  for  e is  proposed  on  solely  dimensional  grounds, 
noting  that  a time  scale  /_1  is  relevant  for  many  dy- 
namical regimes.  Application  of  these  scaling  relations 
in  (10)  yields  the  following  nondimensional  equation, 


. L2  Vh  Ed 
-h+wTz 


(ii) 


Observations  of  turbulence  in  the 
Abyssal  Brazil  Basin 

In  the  period  between  January  1996  and  April  1997, 
two  microstructure  surveys  were  made  in  the  Brazil 
Basin.  The  1996  survey  consisted  of  a basin-scale  sur- 
vey and  the  initiation  of  a tracer  dispersion  experiment 
at  a site  near  the  Mid  Atlantic  Ridge  (MAR).  During 
the  1997  survey,  a detailed  survey  was  conducted  over  a 
5 x 10s  km2  region  of  rough  topography  near  the  tracer 
release  site.  The  primary  focus  of  the  1996  survey  was 
the  spatial  variability  of  mixing  levels,  and  a report  on 
the  observed  variability  is  given  by  Polzin  et  al.  (1997). 
Results  from  the  tracer  release  experiment  are  described 
by  Ledwell  et  al.  (2000). 

The  presentation  here  will  focus  on  the  dissipation 
data  collected  over  the  fracture  zone  (FZ)  topography 
near  the  MAR  (Fig.  2).  All  profiles  consist  of  conven- 
tional hydrographic  variables  (e.g.,  Q,S,U,  V)  in  addi- 
tion to  microstructure,  and  generally  extend  from  the 
surface  mixed  layer  to  within  20  m of  the  bottom.  The 
bathymetric  data  shown  in  Fig.  2 were  derived  from 
satellite  measurements  of  the  marine  gravity  field  Smith 
and  Sandwell  (1997).  The  map  clearly  shows  the  net- 
work of  FZs  leading  east  to  the  MAR.  This  system  of 
FZs  is  characterized  by  a series  of  canyons  bounded 
latitudinally  by  crests  that  rise  up  to  1 km  above  the 
canyon  floors. 

Turbulent  kinetic  energy  dissipation  rates  were  de- 
rived from  observations  of  velocity  microstructure.  Only 
two  components  of  the  strain  tensor  are  measured,  and 
isotropy  is  assumed  to  express  the  dissipation  rate  as 
e = (15/4)-  (<  u2  > + < v2  >).  A summary  of 
the  observed  deep-dissipation  data  was  made  by  col- 
lapsing the  latitudinal  spread  of  the  profiles  (Fig.  3). 
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Figure  2.  Map  of  the  1997  survey  site  with  the  2 arc-minute  resolution  bathymetric  estimates  of  /it  Smith  and 
Sandwell  (1997).  A total  of  129  profiles  of  dissipation  rate  observations  (crosses)  were  collected  in  this  region. 
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Figure  3.  A zonal  section  of  TKE  dissipation  rate  e.  The  individual  profiles  are  shown  as  columns,  and  dissipation 
rate  observations  were  vertically  averaged  into  100  m-bins.  The  shaded  bathymetry  is  representative  of  canyon 
floors,  while  the  shallower  trace  represents  the  level  of  ridge  crests. 
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To  simplify  the  presentation,  only  the  depth  range 
deeper  than  1500  m is  shown.  The  dissipation  data 
were  depth  averaged  into  100-m  intervals,  and  indi- 
vidual profiles  are  shown  as  columns  of  binned  dissi- 
pation. Since  the  latitudinal  extent  of  the  survey  has 
been  collapsed,  no  particular  section  of  bathymetry  is 
entirely  appropriate.  Instead,  two  classes  of  representa- 
tive bathymetry  are  shown:  the  bathymetry  of  FZ  crests 
and  the  bathymetry  of  canyon  floors,  as  sampled  from 
the  latitude  range  21°S  to  23°S.  Within  about  500  m 
of  the  bottom,  there  is  an  enhancement  of  dissipation 
by  more  than  two  orders  of  magnitude  over  values  ob- 
served at  mid  depth.  Enhanced  dissipation  generally 
occurs  between  the  level  of  FZ  crests  and  canyon  floors, 
suggesting  the  largest  mixing  rates  occur  in  canyons. 

An  inverse  estimate  of  circulation 


Surveyed  hydrography  was  used  to  divide  the  flow 
regime  into  a series  of  density  layers.  In  these  layers, 
geostrophic  advection  is  related  to  isobaric  gradients  of 
geopotential  anomaly  (9?)  by 

Vw  . . 

U = z x — -f  Uo,  (13) 

where  u0  is  a reference  level  velocity.  Thus,  the  ad- 
vection of  potential  vorticity  through  a density  layer  of 
thickness  h(x,  y),  as  previously  stated  in  (6),  may  be 
expressed  in  the  form 

lf“o  + (if  - 7)  + £(*».(*,) 

= ur^h  ~ Vt  ("fc'  — f)  ’ 

(M) 

where  (urvT)  are  the  components  of  relative  geostrophic 
velocity  given  by  (u  — uq,v  — vq)  in  (13).  Thus,  the 
problem  of  determining  the  lateral  flow  u is  reduced 
to  determining  the  reference  level  velocity,  i.e.,  the  ver- 
tical integration  constants  (uo,uo)  of  a thermal- wind 
balance. 

In  layers  where  flow  encounters  topography,  the  use 
of  the  geostrophic  balance  is  questionable.  A modified 
momentum  balance,  such  as  that  of  the  bottom  Ekman 
layer,  may  be  called  for  in  regions  where  a density  layer 
comes  within  0(100)  m of  the  bottom.  Additionally, 
the  difficulty  of  defining  the  geopotential  anomaly  along 
pressure  surfaces  interrupted  by  bathymetry  severely 
limits  the  practical  use  of  (13)  and  (14)  in  deep  lay- 
ers. For  these  reasons,  it  is  necessary  to  employ  ad- 
ditional dynamical  constraints  on  the  flow.  This  is 


Figure  4.  Dissipation  profiles  averaged  according  to 
bathymetric  classification.  Bathymetric  classes  were 
assigned  as  crests,  canyons,  or  slopes.  All  profiles 
are  shown  relative  to  a reference  dissipation  of  e = 
1 x 10-10  W kg-1.  The  95%  confidence  interval  is  shown 
for  each  100-m  average. 

achieved  through  the  use  of  integrated  advective  bud- 
gets for  mass, 

1 1 dydz  u\xx+1  + fjdxdzv\yy+A 

(15) 

+ ffdxdyw*  |*“  = 0 

and  for  potential  temperature 


f f dydz  uQT^ 

+ 

f J dxdz  v&\y+A 

+ 

J Jdxdy  w.&'Q 

f f dydz  kQx\1+7 

+ 

f f dxdz  nQy\l+A 

+ 

f jdxdy  kpQz 

(16) 

The  integrated  expressions  (15)  and  (16)  are  bounded 
in  control  volumes  of  meridional  and  zonal  extent  (7,  A), 
and  vertically  bounded  in  a density  layer  by  the  surfaces 
zu(x,y)  and  zt  (x,y). 

Relations  for  the  diapycnal  diabatic  terms  are  ex- 
pressed using  the  dissipation  rate  with  the  use  of  (8)  and 
(9).  For  use  in  the  inverse  model,  the  mixing  efficiency 
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Figure  5.  Layer-averaged  flow  for  the  density  layer  28.16  < an  < 28.20.  The  gridded  surfaces  show  the  an  = 28.16 
(upper)  and  <rn  = 28.20  (lower)  isopycnals.  Bottom  bathymetry  was  subsampled  onto  a 0.25°  x 0.25°  Mercator 
grid.  The  largest  vectors  are  (5-6)  mm  s-1. 


parameter  T is  treated  as  constant,  To  = 0.20  ± 0.04. 
This  allows  us  to  specify  the  diapycnal  diffusivity  for 
turbulence  as 


and  the  diapycnal  advection  as 


(17) 


w * 


To  a? 
N2  dz' 


(18) 


We  emphasize  that  (17)  and  (18)  are  approximations, 
where  a constant  mixing  efficiency  has  been  assumed. 
While  oceanic  observations  of  turbulence  from  well 
stratified  regions  suggest  T0  = 0.2  (Mourn,  1996;  St. 
Laurent  and  Schmitt,  1999),  observations  from  fjords 
suggest  T = 0.05  ( Stigebrandt  and  Aure,  1989).  Ad- 
ditionally, laboratory  measurements  of  turbulence  sug- 
gest 0.05  < T < 0.25  characterizes  the  variation  of  mix- 
ing efficiency  over  a wide  range  Reynolds,  Richardson, 
and  Froude  numbers  ( Ivey  and  Imberger,  1991;  Huq  and 
Britter,  1995).  In  our  inverse  model  estimates,  we  ac- 
cept that  using  To  = 0.2  ± 0.04  may  lead  to  a specifi- 
cation of  kp  that  is  biased  high.  The  use  of  a constant 
mixing  efficiency  in  (18)  is  more  problematic,  as  varia- 
tions in  T with  depth  will  not  only  influence  the  mag- 
nitude of  the  w * estimate,  but  also  the  sign  (w*  > 0 


for  upwelling,  w*  < 0 for  downwelling).  For  this  rea- 
son, to,  is  regarded  as  an  unknown  in  the  inverse  model, 
and  (18)  serves  as  a “constraint”  in  the  inversion.  In 
this  manner,  (18)  provides  an  a priori  estimate,  while 
full  inversion  of  (14),  (15),  and  (16)  determines  the  best 
estimate  of  wt- 

In  expressions  (17)  and  (18),  e is  meant  to  denote 
some  suitably  averaged  function  of  the  dissipation  data. 
Averaging  over  multiple  profiles  is  necessary  to  achieve 
statistical  stability  in  the  mean  dissipation  estimates. 
Since  our  inverse  model  utilizes  the  steady-state  forms 
of  the  advective  budgets,  we  seek  an  averaging  proce- 
dure that  yields  an  estimate  of  the  time-mean  dissipa- 
tion rate.  Thus,  the  dissipation  data  were  averaged  in 
both  space  and  time.  The  time  averaging  was  weighted 
by  a function  of  the  squared  barotropic  tidal-current 
speed  to  remove  biases  associated  with  a spring-neap 
modulation.  Furthermore,  dissipation  data  were  spa- 
tially averaged  according  to  a bathymetric  classification 
scheme.  Dissipation  profiles  above  FZ  canyons,  crests, 
and  the  sloping  topography  between,  were  treated  sep- 
arately in  the  averaging.  The  profile  data,  classified  in 
this  manner,  are  shown  in  Fig.  4.  Each  ensemble  profile 
results  from  data  that  has  been  vertically  averaged  into 
100-m  bins.  Dissipation  decreases  with  height  above 
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Figure  6.  Diapycnal  advection  estimates  for  the  an  = 28.16  (upper  panel)  and  an  = 28.20  (lower  panel)  isopycnals. 
The  lateral  flow  for  this  density  layer  (from  Fig.  5)  is  shown  in  the  middle  panel.  Contoured  estimates  of  diapycnal 
advection  have  been  averaged  over  1°  x 1°  Mercator  cells  which  include  regions  blocked  by  topography. 
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bottom  in  each  of  the  3 classes,  and  all  profiles  reach 
background  levels  of  dissipation,  e = 3 x 10~10  W kg-1, 
at  heights  greater  than  1000  m.  The  slope  profile 
is  notably  different  from  the  other  two,  in  both  the 
magnitude  of  the  maximum  dissipation  and  the  scale 
of  decay.  The  slope  profile  has  a bottom  value  of 
e £ (3  — 9)  x 10 “9  W kg-1  and  decays  with  an  e-fold 
scale  of  (150  ± 50)  m.  The  crest  and  canyon  profiles 
peak  at  e = (2  — 5)  x 10~9  W kg-1  and  decay  over  a 
larger  scale,  (500  ± 100)  m. 

A matrix  model  for  a 6 density  layer  system  con- 
sisting of  about  7000  equations  for  3200  unknowns  was 
formed  using  (14),  (15),  (16),  (17),  and  (18).  The  de- 
tails of  the  linear  inversion  are  given  by  St.  Laurent 
et  al.  (2001).  In  a series  of  initial  calculations,  lateral 
diffusivity  k was  treated  as  constant  along  each  den- 
sity surface  used  in  the  inversion.  Using  this  approach, 
the  estimates  of  k were  inconsistent  with  the  a priori 
estimates  of  k ~ 100  m2  s_1  derived  from  float  obser- 
vations ( Hogg  and  Owens,  1999).  Spatial  variations  of 
k on  isopycnals  intersecting  topography  may  account 
for  this  discrepancy.  Additionally,  diagnostic  scaling 
with  k ~ 100  m2  s-1  indicated  that  a primary  thermal- 
advection  balance  (16)  occurs  between  the  terms  in- 
volving advection  (both  lateral  and  vertical)  and  the 
vertical  diffusion  term.  Scale  estimates  suggest  that 
the  magnitude  of  the  lateral  diffusive  flux  terms  in  (16) 
are  less  than  10%  of  the  diapycnal  diffusive  flux  of  po- 
tential temperature.  For  these  reasons,  terms  involving 
the  lateral  diffusivity  were  dropped  from  the  thermal 
advection  equations,  leaving  the  three  dimensional  flow 
field  (u,  w„)  as  the  model  unknowns. 

Flow  in  and  above  abyssal  canyons 

Here,  we  present  the  estimated  circulation  in  and 
above  the  FZ  canyons.  A more  complete  discussion  of 
the  circulation  at  all  depths  in  this  region  of  the  Brazil 
Basin  is  given  by  St.  Laurent  et  al.  (2001).  In  the  layers 
deeper  than  3000  m,  bottom  bathymetry  has  a clear  in- 
fluence on  the  flow  dynamics.  In  particular,  proximity 
to  the  bottom  is  the  controlling  factor  for  the  magni- 
tude of  the  diabatic  forcing.  The  vertical  diffusivity 
characterizing  these  layers  is  kp  ~ 1 x 10-4  m2  s-1 
or  greater,  and  this  corresponds  to  an  increased  im- 
portance of  the  diabatic  terms  in  the  governing  equa- 
tions. Flow  was  estimated  in  the  neutral  density  layer 
28.16  < <jn  < 28.20  (Fig.  5).  The  upper  isopycnal  for 
this  layer  (ern  = 28.16)  rests  just  above  the  level  of  the 
FZ  crests  at  16°W.  Near  the  MAR,  there  is  westward 
flow  out  of  canyons  with  average  current  magnitudes  of 


(2-4)  mm  s-1.  Further  west,  the  density  layer’s  average 
height  above  bottom  increases,  and  flow  at  500  m above 
the  level  of  FZ  crests  is  southward  with  average  current 
magnitudes  of  (5-6)  mm  s-1. 

Diapycnal  advection  estimates  are  presented  in  Fig.  6. 
These  contoured  fields  of  w * have  been  highly  smoothed. 
Where  FZ  crest  topography  intersects  each  density  sur- 
face, there  is  no  flow,  and  the  contoured  maps  show 
estimates  of  w,  that  have  been  averaged  in  1°  x 1°  cells 
which  often  include  regions  blocked  by  topography.  The 
crn  = 28.16  and  an  = 28.20  isopycnals  are  characterized 
by  diapycnal  downwelling  to  the  west  and  diapycnal  up- 
welling  to  the  east.  Estimates  of  upward  advection  in 
canyons  exceed  w,  = 30  m yr_1,  and  are  as  large  as 
ie*  = 100  m yr-1  in  some  local  regions.  Uncertainty 
estimates  for  the  diapycnal  flow  are  typically  less  than 
±40%  of  the  contoured  values. 

To  produce  a clearer  summary  of  the  deep  circulation 
in  and  above  canyons,  a zonal  stream  function  (ip)  was 
defined  by  latitudinally  integrating  the  velocity  so  that 
uSy  = —ip z and  wlSy  = ipx,  where  the  latitudinal  limits 
of  integration  were  taken  over  the  inverse  domain  (19 °S 
- 25°S), 

u(x,z)  = J dyu(x,y,z), 

_ V (19) 

w.(x,z)  = J dy  w,  (x,  y,  z). 

An  objective  analysis  technique  was  used  to  interpo- 
late the  inverse  solution  and  compute  the  stream  func- 
tion. The  resulting  smoothed  stream-function  is  shown 
in  Fig.  7 relative  to  representative  bathymetry  and  the 
density  field.  These  estimates  suggest  eastward  flow 
in  canyons.  Diapycnal  upwelling  occurs  at  sites  where 
isopycnals  in-crop  along  canyon  floors,  and  this  up- 
welling  feeds  the  westward  flow  at  the  level  of  the  FZ 
crests. 

Deep  upwelling  and  mixing  efficiency 

As  presented  above,  strong  diapycnal  upwelling  is  re- 
quired to  close  the  heat  and  mass  budgets  in  regions 
where  density  surfaces  intersect  canyon  floors.  Figure  8 
presents  a comparison  of  the  a priori  estimates  of  w* 
to  the  inverse  solutions  for  w,  on  the  an  = 28-20  neu- 
tral surface.  Inverse  estimates  are  generally  consistent 
with  the  a priori  estimates  of  w,  at  sites  where  the 
an  = 28.20  surface  is  more  than  300  m above  the  bot- 
tom. However,  the  inverse  estimates  suggest  diapycnal 
upwelling  is  favored  at  sites  where  the  an  = 28.20  sur- 
face is  within  300  m of  the  bottom. 
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Figure  7.  The  meridionally  integrated  stream  function  as  estimated  through  objective  analysis  of  the  inverse 
model  solution.  Dimensional  values  of  the  stream  function  are  contoured,  and  uncertainties  are  roughly  ±40%. 
Representative  bathymetry  of  the  canyon  floors  (shaded)  and  crests  (gray  line)  is  shown,  as  are  the  an  — 28.10,  28.16 
and  28.20  isopycnals. 


Diapycnal  upwelling  must  be  occurring  in  regions 
where  turbulence  supports  a convergent  buoyancy  flux: 
w*N2  = —dJb/dz  > 0.  This  condition  is  satisfied  by 
turbulence  occurring  close  to  the  bottom,  since  the 
buoyancy  flux  must  decrease  with  depth  to  meet  a no- 
flux condition  at  the  sea  floor.  Within  the  context  of 
(7),  the  buoyancy  flux  is  related  to  the  dissipation  rate 
by  a mixing  efficiency  parameter,  Jb  = — Te.  We  believe 
that  a reduction  in  mixing  efficiency  near  the  bottom 
allows  the  no-flux  condition  to  be  met  at  the  sea  floor. 
This  reduction  may  occur  in  a 1 to  10  m thick  bound- 
ary layer  (i.e.,  the  “log  layer”),  which  was  unresolved 
by  our  measurements.  However,  variations  in  mixing 
efficiency  may  extend  to  greater  heights  in  the  canyons. 
The  physical  mechanisms  controlling  the  efficiency  of 
mixing  near  topography  are  not  well  understood,  but 
turbulent  mixing  in  regions  enclosed  by  topography  may 
favor  mixing  with  reduced  efficiency.  Stigebrandt  and 
Aure  (1989)  found  that  mixing  in  fjords  occurred  at  0.05 
efficiency.  Given  our  observations  of  dissipation  and 
stratification  along  canyon  slopes,  a decrease  in  mixing 
efficiency  from  0.2  to  0.05  in  the  bottom  most  100  m 
above  FZ  slopes  would  account  for  diapycnal  upwelling 
as  large  as  w*  = +300  m yr-1. 


Discussion 

We  have  described  buoyancy  forcing  and  circulation 
occurring  in  a region  of  the  abyssal  ocean.  Our  cal- 
culations suggest  that  enhanced  levels  of  turbulence 
above  rough  bathymetry  lead  to  enhanced  levels  of  di- 
abatic  flow.  Divergence  of  the  diabatic  flow  acts  as  a 
mechanism  of  buoyancy  forcing  for  the  lateral  circula- 
tion though  the  vortex  stretching  term  (/ /h)(w*(zu)  — 
w*(zi ))  in  the  vorticity  budget.  The  significance  of  di- 
abatic vortex  stretching  as  a forcing  mechanism  must 
be  assessed  relative  to  other  forcing  agents  for  the 
flow.  From  the  inverse  estimates,  we  may  compare  the 
strength  of  buoyancy  forcing  to  the  adiabatic  mecha- 
nism of  vortex  stretching,  represented  by  ( f/h)u  ■ Vh 
in  the  vorticity  budget. 

Above  the  level  of  FZ  crests,  adiabatic  vortex  stretch- 
ing is  cyclonic,  while  diabatic  vortex  stretching  is  anti- 
cyclonic.  In  the  region  where  mixing  rates  are  largest, 
these  two  mechanisms  are  comparable  in  magnitude  and 
the  vortex  stretching  terms  in  (6)  nearly  cancel.  Since 
/3v  = 0,  there  is  little  meridional  flow  just  above  the 
level  of  FZ  crests  (Fig.  5).  The  divergence  of  diapycnal 
mass  flux  is  clearly  a significant  forcing  agent  on  the 
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Figure  8.  Comparison  of  the  a priori  estimates  of  wm  with  the  inverse  estimates  for  the  an  = 28.20  neutral  density 
surface.  Estimates  at  sites  where  the  <r„  = 28.20  surface  is  more  than  300  m above  the  bottom  are  plotted  using 
triangles,  and  estimates  from  sites  where  an  = 28.20  is  less  than  300  m above  the  bottom  are  plotted  as  points. 
The  curve  along  which  the  inverse  and  a priori  estimates  are  equal  is  shown.  Standard  errors  are  shown  for  all 
estimates. 


lateral  circulation  in  this  region.  Further  west  of  the 
Mid  Atlantic  Ridge  where  mixing  rates  are  weaker,  the 
adiabatic  stretching  exceeds  the  diabatic  stretching,  re- 
sulting in  the  southward  flow  where  the  largest  lateral 
velocities  occur.  In  addition  to  the  stretching  mecha- 
nisms, friction  imposed  by  topography  may  also  be  a 
significant  forcing  mechanism  near  the  bottom,  though 
we  have  not  attempted  to  assess  this  contribution. 

The  diapycnal  advection  estimates  presented  here 
are  clearly  sensitive  to  assumptions  made  about  tur- 
bulent mixing  efficiency.  As  indicated  by  (9),  depth 
variations  in  mixing  efficiency  will  influence  both  the 
magnitude  and  direction  of  the  diabatic  mass  flux.  We 
note  that  most  data  are  consistent  with  a mixing  ef- 
ficiency Rf  near  15%  (equivalently,  T = 0.20).  For 
example,  diffusivity  estimates  from  (17),  which  assume 
T = 0.2  ± 0.04  agree  with  the  turbulent  diffusivity  es- 
timated from  the  tracer  dispersion  experiment  ( Ledwell 
et  al,  2000).  Moreover,  St.  Laurent  et  al.  (2001)  find 
general  agreement  between  a priori  diapycnal  advection 
estimates  estimated  with  F = 0.20  in  (18)  and  diapyc- 
nal advection  estimates  deduced  from  inversion  of  the 
steady-state  budgets.  Variations  in  mixing  efficiency 


are  only  inferred  from  estimates  deep  in  FZ  canyons,  in 
regions  where  isopycnal  slopes  increase  as  density  lay- 
ers in-crop  along  the  bottom.  Our  evidence  for  reduced 
mixing  efficiencies  is  indirect,  but  the  diapycnal  ad- 
vection estimates  are  consistent  with  a reduction  from 
T = 0.20  to  T = 0.05  over  a depth  of  100  m above  the 
sloping  walls  of  FZ  canyons. 

While  the  study  strictly  considered  data  from  a small 
geographic  region,  the  bathymetric  features  present  are 
ubiquitous  over  the  basins  of  the  global  ocean,  as  are 
the  tides.  If,  as  suggested  by  Munk  and  Wunsch  (1998), 
the  tides  are  the  dominant  source  of  mechanical  energy 
for  abyssal  mixing,  the  pattern  of  enhanced  turbulence 
observed  over  Brazil  Basin  fracture  zones  may  be  repre- 
sentative of  turbulence  occurring  in  other  regions.  Dia- 
batic forcing  above  rough  topography  may  be  a primary 
forcing  mechanism  for  the  abyssal  circulation. 
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